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O. Abstract 

In this paper a spin-1 spin-glass model under the presence of a uni- 
^ ' form crystal field is investigated. It is shown that the model presents 

■^ . both continuous and first-order phase transition separated by a tricriti- 

cal point. The phase diagram is obtained within the replica-symmetric 
ff^ ■ solution and exhibits reentrance phenomena at low temperatures. Pos- 

1/-^ . sibly it is the simplest model which can describe inverse freezing phe- 

^D ' nomena. 

o 



Tricritical behavior and reentrance phenomena associated with first-order 
^ I transitions in disordered systems have been the subject of recent studies. For 

instance, investigations of inverse freezing phenomena, where the ordered 
phase is more entropic than the disordered one, have been conducted by 
several authors [H [21 [3l HI O [6]. Such nonusual behavior has also been 
observed in fermionic disordered systems (see [7] and references therein). 
More recently, a spin model for strain glass has been used to account for the 
effects of disorder in ferroelastics which can display reentrance phenomena 
[8]. The models involved in those investigations are closely related to the 



Ghatak-Sherrington (GS) model [H], which has been intensively studied since 
its introduction [U UHl HIl [HI [121 E] • 

The Ghatak-Sherrington model [9] is a generalization of the widely known 
infinite-range Sherrington-Kirkpatrick (SK) model for a spin glass (SG) ^15] 
with arbitrary spin S > 1/2 and inclusion of a uniform crystal field. For 
integer spin 5* the GS model displays both first-order and continuous tran- 
sitions. The crystal field (D)-temperature (T) phase diagram for 5* = 1 has 
a continous transition line which meets a first-order transitions (FOT) line 
at a tricritical point. These lines separates the paramagnetic from the spin- 
glass phase. One of the most interesting feature of the S = 1 SG model is 
the appearance of reentrance effects occuring at low temperatures which has 
been recently associated with inverse freezing phenomena [H [21 El IH El E] • 

Although it is widely accepted that the correct mean-field solution for 
infinite-range spin-glass models is given by the Parisi Ansatz [HI [T71 [18] , the 
replica-symmetric (RS) solution gives an initial clue to general topologies of 
phase diagrams. For the GS model the RS solution has revelead a number 
of difficulties in the low temperature region where first-order transtions take 
place [12j. Inside this region we can find up to three distincts paramagnetic 
solutions. However, only one of these solutions is stable with respect to 
replica-symmetry breaking (RSB) fiuctuations. On the other hand, we may 
find up to four distincts spin-glass solutions in the same region which, unlike 
the paramagnetic case, are all unstable with respect to RSB. Therefore, even 
within the RS solution treatment we are faced with the problem of choosing 
the most adequate SG solution in order to determine the FOT line. We 
believe that this fact is at the origin of some controversies with respect to 
the location of the FOT line predicted by the RS solution [9l [TOl [11] . Recently 
we have seen some progress with respect to the correct location of the FOT 
line for the GS model P^ [Tl] . There remain, however, some points which 
need to be clarified about the low temperature behavior of the GS model. 

Another route of investigation on the mean-field behavior of spin-glasses 
is the Thouless-Anderson-Palmer (TAP) approach which avoids the replica 
method [121 . From the TAP equations one can obtain a simpler set of equa- 
tions by excluding the so-called Onsager reaction-field term so that one gets 
the naive mean-field equations. Bray, Sompolinsky and Yu (BSY) [20j intro- 
duced an exactly soluble infinite-range SG model which is also a generaliza- 
tion of the SK model. For the BSY model one can obtain exactly the naive 
mean-field equations. These authors also studied their model by means of 
the replica method and noted some interesting features. For instance, the 



RS solution describes a SG phase whose entropy is always non-negative for 
finite temperatures and vanishes at zero temperature, in contrast to what is 
observed in the RS solution to the SK model. In spite of this, the RS solution 
to the BSY model is unstable in the whole SG region. Due to lack of the 
Onsager term the BSY predicts a critical temperature that is twice the result 
obtained for the SK model. Notwithstanding this, the low temperature be- 
havior of both BSY and SK models are qualitatively the same. In particular 
they have the same zero-tempareture properties. 

In order to gain further understanding about tricritical behavior and reen- 
trant effects in spin-glass systems we consider a BSY version of the GS model. 
Our analysis is based on the replica approach, since we want to compare our 
findings with the known results obtained for the GS model in previous stud- 
ies. 

The model consists of a set of m classical spin-1 variables Sia (a = 
1, . . . , m) located at each site i = 1, . . . ,N. The Hamiltonian is given by 

^ m N m 

^ = -^ E E J^^^-^^^ + ^ E E ^-' (1) 

(j,jf) a, 6=1 j=l a=l 

where Sia = ±1,0; the {i,j) sum is over all distinct pairs of sites; the ex- 
change interactions Jij are quenched random variables with the Gaussian 
distribution 

and D represents the effect of a uniform crystal field anisotropy term. Several 
known cases are recovered in special limits: (i) m = 1 recovers the Ghatak- 
Sherrington model; (ii) for D — )> — oo the BSY model is re-obtained, including 
the SK model for m = 1. 

Since we are interested in the m — )■ oo limit, the quenched free energy 
per spin is given by 

-(3f= hm hm -^ (In Z) j , (3) 

where ()j denotes the configuration average over the disorder. In order to 
proceed further the use of the replica method is introduced through the 
identity InZ = lim„_^o(^"'^l)/'^- Therefore, the averaged free energy density 
may be expressed as 



where 



/3/ = lim lim min (pmniiqaiB}) , (4) 



1 " 

i({9"/3}) = -^ni^iPJ? XI ^^/3 ~ lnTr{s.}exp(7/^„), (5) 

a,/3=l 



and 



n m 



y-mn = lipj)' E E ^"z^^"^^'' - ^^ J2(s:r. (6) 



2 

a,/3=l a, 6=1 a,o 



The condition for ^mnd^a^}) to be an extremum with respect to qa/s 
yields 

^°/3 = ;;^ E (^aSb)mn (7) 

a,fe=l 

where {S2S^)mn indicates the thermal average with respect to the replica 
hamiltonian (E]). 

The replica-symmetric solution can be considered now. As in the BSY 
model there is a non-trivial diagonal qa^ and the requirement of a finite 
susceptibility leads us to consider the RS Ansatz in the form: 

qap = q (a 7^ 13), qaa = q + x/m. (8) 

Substituting ([8]) into (jl]) the free energy per spin is obtained 

(3f = ^{PJfxq-(-li(3Jfxm'ix) + lnzix)\ , (9) 



where 



and 



zix] 



1 + 2e-^^ cosh{pj^x + (3^J^xm{x)), (10) 



^ 2 sinh(/3 jygx + /3V^x"-^(3^)) ,^^x 

"^^^'' - g/3D + 2 cosh(/3 J^x + P^J^xm{x)) ' ^ ' 



The notation has been simphfied by introducing 



^27r 
The equihbrium equations df/dq = = df /dx yield 

q={m\x)),, (12) 

I3x = ^{xm{x))^. (13) 

An integration by parts allows to re-write ( 1T3|l as 

^-, _/ p{x)-m^{x) \ 

^ \l-(/3J)2^(p(a;)-m2(x))/,' ^ ' 

where 

2 cosh(/3 J^x + (3^J^xmix)) 
P^^> ~ eP^ + 2 cosh(/3 jyga; + P^J^xH^)) ' 

The set of Eqs. iQ, (TT2|) and ( [T3|) determine the phase diagram which 
presents paramagnetic and spin-glass phases. It should be mentioned that 
for numerical purpose Eq. (fT4l) is more appropriate than Eq. (TT3!) . 

The paramagnetic phase is described by g = and 



The simple form of p in the last equation shows that it is a single val- 
ued function both in terms of inverse temperature /3 and anisotropy crystal 
field D. It is important to recall that in the Ghatak-Sherrington model the 
behavior of the corresponding term is more subtle [12]. As a matter of fact 
in the GS case p represents a true order parameter and can display up to 
three distinct solutions. Thus we have to find an additional criterion in or- 
der to determine the thermodynamically stable paramagnetic solution. In 
the present case p can be regarded as a mere density and a simple analysis 
shows that the paramagnetic solution is physically acceptable as long as the 
condition 



'''-'' < 1 (17) 



ePD + 2 
holds. 

Let us set J = 1 and consider the D — T phase diagram, where T = 
1//3. From Eq. (ITTj) the paramagnetic is stable at high temperatures and is 
bordered by the line 



D = Tln 



2(2 -T) 



(18) 



At low temperatures there is a spin-glass phase with g > 0. Expanding 
Eq. flT2l) for small q one finds 

q = aq + bq"^ + 0{q^) (19) 

where 



"^ (20) 



and 



T"^ -Xp 



l(l-3p)^ . 



^ 9 .3 Ta^ (21) 

Thus the spin-glass phase exists as long as a < 1 and b > 0. From the 
above expansion one finds a tri critical point given by a = 1 and 6 = 0. 
Therefore the phase diagram consists of a continuous transition line given by 
Eq. ( !T8|) as long as T > 2/3. For < T < 2/3 there is a region of coexisting 
paramagnetic and spin-glass solutions and one has a first-order transition 
which can be numerically determined by equating the corresponding free- 
energy densities of these solutions. The continuous and first-order transition 
lines meet at the tricritical point given by 

T = 2/3, D = 4 In 2/3 = .924169..., (22) 

which should be compared with the correponding results for the tricritical 
point found for the GS model: T = 1/3, D = l/2 + 2(ln2)/3 = 0.962098 . . . 

i- 

At zero temperature the first-order transition can also be easily found. 

First one notes that \imT-^o /^X = X is finite. From Eq. (IT^ one finds 

6 



oo 

.2/ 



q = 2 I e-"'/2^ = eTic{x*/V2), (23) 

27r 



where 

and erfc is the usual complementary error function. The above expression is 
valid for x* > 0. In this regime one also finds 

X^(A) e--/.. p6) 

and, for the spin-glass free-energy density /o, 



/o = Dg - (^^J e-^' /2. (26) 

For negative values of x* we obtain, at the absolute zero of temperature, 
q = 1, X = (2/7r) and / = — (7r/2) / . Since from Eq. (|9]) the free energy 
of the paramagnetic solution at T = is zero for D > 0, the first-order 
transitions at T = is found by imposing, /o = from which follows, after 
some simplifications with the help of Eqs. ( 123|) and ( l25l) 



X* eMx*/V2) = ^e""* ^^ (27) 

V 27r 

The above equation is exactly the same obtained previously for the replica 
symmetric solution [12] to the GS model. The numerical solution to this 
equation is x* = 0.612003 . . . , from which results, jointly with Eqs. (123|) and 
i^^, q = 0.540535 . . . and x = 0.899003 . . . , respectively. These numerical 
results allow us to determine the location of the first-order transition at 
T = 0: 

Do = D{T = 0) = 0.899033 .... (28) 

Since the present model as well as the GS model have the same ground 
state this is an expected result of general validity. Thus, as in the GS model 
a stable RSB solution should give a slightly lower value for D^-, and so the 
reentrance effect must be enhanced. In fact, a previous numerical study of 
the naive mean-field equations for T = showed that Dq ^ 0.86 



For < T < 2/3 the first-order transition line can be obtained by nu- 
merically solving the equilibrium equation and requiring that the spin-glass 
and paramagnetic solution have the same free energy. The resulting phase 
diagram is depicted in Fig. 1. As in the GS model case there is a reentrance 
to the paramagnetic phase at low temperatures. However, one notices that 
as T — > there is no new transition to the SG phase in the vicinity of Dq 
as was found in the replica solution to the Ghatak-Sherrington model [T^ . 
This is a direct consequence of the vanishing of the spin-glass entropy as can 
be verified by an analysis of the Clausius-Clayperon equation along the FOT 
line. 

In spite of the vanishing of the spin-glass entropy at T = for any 
value oi D, a stability analysis along the lines pioneered by de Almeida and 
Thouless [2T] shows that the replica-symmetric spin-glass solution is always 
unstable. Again, the numerical analysis of the stability conditions becomes 
easier than in the corresponding Ghatak-Sherrington case since we did not 
find any evidence of complex eigenvalues for the replica stability matrix. The 
instability is signaled by the non-positiveness of the replicon eigenvalue 

XR = T'-{{p{x)-m\x))'),. (29) 

Therefore, a correct description of the spin-glass phase requires a complete 
solution to the corresponding Parisi's equations as have been done for the GS 
model [131 [H]. We can anticipate that in such replica-symmetry breaking 
treatment the main features of the phase diagram would not be modified. 
The continuous transition line as well as the location of the tricritical point 
will not be changed, but a slight modification in the location of the FOT line 
is expected in a similar way to what happens in the GS model, increasing 
teh reentrance effect for T ^ 0. 

In conclusion, a naive version of the Ghatak-Sherrington model was inves- 
tigated by the replica approach. The phase diagram was determined within 
the replica-symmetric Ansatz. Both analytical and numerical results shows 
that the replica-symmetric solution to the present model is simpler to anal- 
yse than the corresponding Ghatak-Sherrington model. In spite of this, the 
phase diagram for both models share several common features. At high tem- 
peratures there is continuous transition line from the paramagnetic to the 
spin-glass phase, while at low temperatures there is a line of first-order tran- 
sitions. These two lines meet at a tricritical point. The present model also 
exhibits a reentrance from the spin-glass to the paramagnetic phase. We 



believe that the present model could be useful to further investigations on 
inverse freezing phenomena. 

The author thanks to Prof. Silvio Salinas for useful comments. 
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Figure 1 Phase diagram obtained within the rephca-symmetric approxima- 
tion. The full AT curve is the line of continuous transitions and the broken curve 
BT is the line of first-order transitions. These two curves meet at the tricritical 
point T. 
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